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O ' Abstract 

^ . We propose a simple phenomenological parameterization of quint- 

essence with a time- varying equation of state. In particular, it ac- 
counts for the possibility of early dark energy. The quintessence po- 
tential can be reconstructed in terms of the present fraction in dark 
energy, the present equation of state and the amount of early dark 
energy. 

> 
ON 

Distinguishing quintessence - a time varying dark energy component - 
■ from a cosmological constant is a major quest of present observational cos- 

mology. The first particle physics models jH |2] are based on a scalar field 
rolling down a potential. Typically, a cosmic attractor solution |UI2] renders 
the evolution of the scalar field independent of the initial conditions ( after a 
certain transition time in early cosmology). For a given scalar potential the 
quintessence cosmology is as predictive as a cosmological constant scenario. 
However, there are numerous models with different potentials. For cosmolog- 
ical tests of their general features one would like a simple parameterization 
of the cosmological dynamics, say the Hubble parameter as a function of 
redshift, H(z), that can later be translated into statements about the form 
of the quintessence potential. 

A first obviously important parameter is the present fraction in homoge- 
neous dark energy pP £1® = Vth{z = 0). For observations at low z the next 
important parameter is the first derivative of flh(z) or the present equation 
of state wq = Wh(z = 0) [3 . These two quantities are directly related, cf. eq. 
([3*]) below. Continuing an expansion in z is not very meaningful, however, 
if one wants to cover the physics at high z, like last scattering. For exam- 
ple, the CMB-anisotropies are very sensitive to the fraction of dark energy 
at last scattering jlj = Vt h {z = 1100) and structure formation depends 
crucially on the weighted mean of Qh during structure formation IE] • 
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These two important parameters are not reasonably described by a Taylor 
expansion of a h (z) or Wh(z) around z — 0. 

The long term goal is certainly to gather as much information as possi- 
ble about the whole function Qh(z) [Zj. One possible strategy leaves flh(z) 
essentially unconstrained and simply tries to find a function which improves 
substantially the fit to the available data as compared to a cosmological con- 
stant. A complementary strategy describes £lh{ z ) [HE] or Wh(z) HE EH HI| 
in terms of a few parameters. In particular, the parameterization of Wh by 
Fermi-Dirac functions [10] can cover many interesting proposed quintesssence 
models. We suggest here a concentration on Qh and on the absolute minimum 
of parameters: we propose to test a three parameter family of quintessence 
models characterized (for a flat universe) by a M = 1 — , w and a 
new bending parameter b or, equivalently, the fraction of early dark en- 
ergy, Q e . The criteria for the parameterization are chosen such that at least 
for a certain parameter range the model is consistent with a scalar field 
rolling down a potential. (This is not the case, for example, for an ansatz 
w(z) = w + w'z except for w' = 0.) It also should cover the interesting 
case of early quintessence. The two parameters 1 + Wq and b measure the 
deviation from the case of a cosmological constant. At least in principle they 
can be determined from every cosmological observation separately since our 
parameterization covers the whole available redshift range. If, at the end, 
values of 1 + Wo and b different from zero are preferred and their best values 
differ between observations covering different redshift ranges, it will still be 
time to enlarge our parameterization. 

Three parameters 

Observations have most direct access to the redshift dependence of the 
Hubble parameter, H(z). The dark energy quantity that is most directly 
related to H(z) is the fraction in dark energy as a function of redshift, Qh(z). 
We aim here at a useful parameterization of f^(z) in terms of three (or four) 
parameters. Two parameters are chosen as Qm — 1 — ^i '' and wq. In order 
to proceed further we use a simple relation between Qh(z) and Wh(z) which 
is valid if the energy density not contained in dark energy can be accounted 
for by pressureless matter [TT| 
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This suggests to parameterize the function 
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which obeys (in absence of radiation) 



dy 

We propose 



9R(y) Q * v /q\ 
"= 3tw ft (y). (3) 



*(,)=* + £| (4) 



where R is directly related to VL M by 



*-(Hsr)- (5) 

The new "bending parameter" b = 1/yb = l/ln(l + Zb) characterizes the 
redshift where an approximately constant equation of state turns over to a 
different behavior. Nonzero b signals the breakdown of the linear expansion 
for R(y). Within the simple parameterization (j3J a positive b is directly 
related to the presence of early dark energy 

n n i \ exp(fi + 3w /b) 

Q e = n h (y ^00 = — — — — . 6 

1 + exp(i? + 3w /b) 

Whereas for supernovae (and other observations at low z) b is perhaps the 
most natural parameter we actually suggest to quote the triplet {VL M , w , Q e ). 
The parameter Q e has a simple physics interpretation. For practical purposes 
within the three parameter parameterization Q e equals and fijf and 
enters therefore very directly the CMB anisotropics or structure formation. 
Inserting eq. (jOJ b is then determined as 



3u>n 



In 



^)+ln(^ 



(7) 



(Nevertheless, for the fitting to the data b is preferred since eq. (j3J) also 
covers the case of negative b. Supernovae observe moderate values of y ~ 1 
and the parameter b may be related to Vth{z = 1) which is perhaps most 
directly probed.) For a cosmological constant one has wo = — l,Q e = or 
wq = — 1, b = whereas a constant equation of state obtains for 6 = 0. If 
one performs a Taylor expansion of Wh{z) [S] 

w h (z) = w + w'z + . . . (8) 

one finds the relation 

w' = -2w b. (9) 
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Neglecting radiation, the function R(y) has a simple interpretation in 
terms of the averaged equation of state 



y 

My) = -( dy'w h (y') = R{y) - Ro . (io) 

y J 3y 

o 

This yields a very convenient formula for the redshift dependence of the 
Hubble parameter according to 

= (i - n M )(i + + q m (i + zf. (ii) 

Therefore the luminosity distance and related cosmological quantities are 
easily expressed in terms of our parameterization of R. With (jH) one finds 

= T+m+Ty (12) 

Our ansatz can also be interpreted as a direct parameterization of the redshift 
dependence of the Hubble parameter. 

We note that our parameterization ensures automatically that flh(z) re- 
mains between and 1 for arbitrary values of R(y). It implies a negative 
equation of state (for w < 0) 

Mv) = jY^yr (13) 

Therefore R(y) and Qh(z) are monotonic functions. For quintessence due to 
the evolution of a scalar field one has 

T - V 

w h = (14) 

h T+ y V ) 

with kinetic and potential energy T and V. For positive T and V the equation 
of state is bounded to the interval — 1 < < 1. (For T > 0, V < one 
finds \wh\ > 1, with sign depending on the sign of T + V.) One may argue 
that the form (|l~3|) is not general enough and, in particular, not suitable for 
the description of a possible change of sign of Wh- In fact, for typical models 
with early dark energy one expects Wh > in the radiation dominated epoch. 
We remark, however, that eqs. ©,© receive corrections ^T] once radiation 
becomes important. (For Q e > one finds for the radiation dominated epoch 
Wh = 1/3.) Our parameterization can therefore remain reasonable as long 
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as Flh( z ) is a monotonic function 2 which is the case for most quintessence 
models. 



Reconstructing the scalar potential 

The scalar potential and kinetic term can be reconstructed ^Tj from R(y). 
As a consequence, our three parameters describe a family of complete models 
for which not only the background evolution but also the fluctuations of the 
cosmon field can be determined. One first constructs V(y) according to 



v= 1 -^ Ph = ^f(i-w h )n h H* (15) 

(with reduced Planck mass M 2 = M 2 /8tt, p cr = 3M 2 H 2 ). For the relation 
between the scalar field ip and y we also use the scalar kinetic energy 

t = —(i + Wh )n h H 2 

One possibility employs a standard kinetic term (k = 1) and integrates the 
relation 

1 d(p 



MTy = ^(l + ^R ( 17 ) 

in order to obtain <p(y), and, with ()15|) . V(f). Perhaps more direct and 
more elegant is a rescaling (p = <p((j)) such that the potential V((f>) takes a 
standard form. The information about the specific quintessence model is 
then contained in the "kinetial" k(<f>) [I2j- In this case <p(y) can be directly 
extracted from eq. (JT3J). Subsequently, eq. (fTBj) can be solved for k. In 
particular, if the potential is monotonic in the scalar field one may use the 
"standard exponential form" 



V = M*ex V ^) (18) 

For this choice one extracts from eqs. (12) (14) of ref. jTI] 

AT 1 = [3(1 + w h )Q h ]- l/2 (3(1 + w h ) - ^_J— ) (19) 
I oy 1 - w h J 

2 Actually, in many models the potential is essentially exponential at early times. In 
this case a calculable change £Q in Clh(z) occurs at the transition from the radiation to 
the matter dominated era. Beyond Sl e this jump does not involve a new free parameter 
and could easily be incorporated into our parameterization. Since we concentrate here on 
the matter dominated epoch we have left out this issue in this note. 
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and from eqs. (|T5Jl (fTHj) 



-£ = -\ n (HL(i- Wh )n h \ (20) 



In conclusion, we can associate to a given parameter set w ,b) a quintessence 
model with a specific kinetical /c(0), or after rescaling to a standard kinetic 
term, a potential V^(<^). 

A few comments illustrate this relation. 

(i) For O;,>0a positive T requires Wh > — 1. 

(ii) For a region with > and > 1 the potential V must be negative 
(cf. eq. (I15|) ). The parameterization (jlj) does not cover this case. 

(iii) The relation (JT5j) holds only up to a sign. For k~ l = the potential 
V(ip) has a stationary point, dV/dip = 0. A change of sign of k~ l 
according to eq. ()19|1 signals that the scalar field moves through an 
extremum of the potential V(ip) at the corresponding value of y. This 
implies an interesting condition for a monotonic V(ip), namely 3 



dy 



^3(l-wi). (21) 



Within the ansatz Q k 1 is positive for large y if b > 0. A monotonic 
potential therefore requires for y = 

or (wo < 0) 

6 < - 3 ' 1 -°>. (23) 
2w„ k ' 

Our parameterization covers a wide class of dynamical behaviors of dark 
energy beyond the quintessence models based on a scalar field. We recall 
that w < 1 or b < are not compatible with scalar quintessence models 
(with T > 0). For w close to —1, as suggested by observation, only a rather 
narrow range of b is compatible with a monotonic potential. As an immediate 
consequence, quintessence models with monotonic potential and leading to 
Wo close to —1 are very hard to distinguish from a cosmological constant 

3 This condition follows also directly from eq. (|13|) . It is equivalent to dlnV/dy ^ 0, 
using d\np h /dy = 3(1 + w h ). 
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by observations at low z. Indeed, w(z) is bounded by the solution of the 
differential equation dwh/dy = 3(1 — w%) (cf. eq. (jUJ)), namely 

_ i- Wo -(i + Wo yy 
Wh[y) ~ i- Wo + (i + w yy- {2A) 

Independently of our specific parameterization a monotonic potential requires 
Wh(z) < w^ p (z). For small (1 + Wq) one needs sufficiently large z before 
w(z) can deviate substantially from —1. Actually, the bound ({21} cannot 
be saturated since this would require that the scalar field is frozen at an 
extremum, in contradiction to Wh > — 1. Combining the bound (|2*3*j) with eq. 
(|13|1 yields a more severe bound within our parameterization 

Whiy) < — q ~ 9 (25) 

' *£2& y Y [1 + 3(1 + w )y] 2 



2w 



where the last expression assumes 1+Wo < 1- In this case, i.e. for b obeying 
the inequality (}23*j) . we find a monotonic fc(0). 

Extended parameterization 

Combining the inequality ()23|) with eq. (J2J) restricts fi e to rather small 
values if 1 + wo <C 1. This is, of course, a consequence of our ansatz (J3J) 
with only three independent parameters. From the particle physics point of 
view there is no contradiction between a sizeable fl e (say a few percent) and 
Wo close to —1. One may therefore want to weaken the strict connection 
between b and Q e by an extended parameterization. For example, this may 
be required if supernovae find a small b while structure formation or the CMB 
prefer a nonvanishing Q e . A still relatively simple possibility is the choice 

nh(y)= (i-n Ar -n e )-i +n - (26) 

where 

(i-^expfe^) 

n«(y = / \ (27) 

+ (l - n*) «p (^t^J 

coincides with Qh{y) and the ansatz (0J) for a = 1. The relation between b 
and f2 e depends now on the new parameter a 



3w 



:28i 
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For a Taylor expansion Wiz) = Wq + w'z one finds 



3(a - 1) Q 2 h - QJ2Q h - 1) 9 , , 

u;' = -26 Wo + ^ >- h e{ h '-wl (29) 

a \l h — \l e 

Another possible shortcoming of the ansatz (jlj (and also the extension 
()26)) ) is the necessarily nonzero value of fl e for b > 0. One may include the 
possibility that Qh does not approach a constant for large y even for b ^ 
and use an alternative extended parameterization 

^) = Ro + — — (30) 

with 

w + bjy{2 + by) 

My) = {1 + by)2 ■ (3i) 

In practice, one may express 7 in terms of the fraction in dark energy at last 
scattering, Qf^ = Vth{z = 1100). Negative 7 implies that Qh approaches zero 
for y — > 00, whereas it formally increases towards one for 7 > 0. Neglecting 
radiation 7 denotes the "asymptotic" value of Wh(z — > 00) = 7. (For this 
particular point we assume b > such that R has no pole for y > 0.) Exten- 
sions of the parameterization like 1)26)1 or (jHH)) may be needed in the long run 
for a detailed comparison of CMB, structure formation and supernovae. For 
observations covering a range y ^ 1 the parameterization (jlj) is presumably 
sufficient. 

Other parameterizations of a time varying equation of state Wh(z) or the 
corresponding flh(z) have been proposed in the literature |U El El El EEl 
ITT] . They often cover only a restricted range of redshift but fail (or are 
unspecified) when applied to the full redshift range accessible to observation, 

say < z ^5 1200. For example, the Taylor expansion w(z) = wo + w'z leads 

to w(z) > 1 for z > z cr = (1 — w )/w' {w' > 0). This only is consistent with 
a very restricted and not very natural class of scalar models, namely those 
where the potential V(<p) was negative in the past for z > z cr and has turned 
positive recently for z < z cr . (For negative w' the positivity of kinetic energy 
{w(z) > 1) is violated in the past - even worse.) On the other hand, the 
proposed parameterizations covering the whole redshift range lead usually to 
a more complicated form of H(z) as compared to eq. ()TT)) . involve more than 
three parameters and often contain a certain degree of arbitrariness. 

It is obvious that every parameterization in terms of a few parameters 
will have its strength for a particular problem but also its shortcomings when 
applied to the whole relevant redshift-range. Nevertheless, concentrating on 
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a simple one would be beneficial for a direct comparison of different obser- 
vations at the present stage. 
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